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Abstract. If the Higgs boson H couples to a singlet scalar S via Xm\H\'^S'^, a strong elec- 
troweak phase transition can be induced through a large potential barrier that exists already 
at zero temperature. In this case properties of the phase transition can be computed ana- 
lytically. We show that electroweak baryogenesis can be achieved using CP violation from a 
dimension-6 operator that couples S to the top-quark mass, suppressed by a new physics scale 
that can be well above 1 TeV. Moreover the singlet is a dark matter candidate whose relic 
density is < 3% of the total dark matter density, but which nevertheless interacts strongly 
enough with nuclei (through Higgs exchange) to be just below the current XENONIOO limits. 
The DM mass is predicted to be in the range 80 — 160 GeV. 
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1 Introduction 

New physics is required for viable electroweak baryogenesis and for dark matter (DM). A 
natural scenario for potentially achieving both simultaneously is scalar dark matter S that 
couples to the Higgs field H through the renormalizable interaction Xjn\H\'^S'^. The case 
where the new scalar is another SU(2)i doublet (the inert doublet model model [1]-[16]) 
has recently been studied in this context [17]-[19], where working examples were found, at 
the expense of some tuning of parameters in order to sufficiently suppress the relic density 
to satisfy direct detection constraints while getting a strong electroweak phase transition 
(EWPT).^ The case of singlet S dark matter [20]- [34] has also been considered with respect 
to its impact on the electroweak phase transition in refs. [35]- [39] (for the effect of singlets 
on the EWPT without requiring them to be dark matter, see [40]-[48]). It was found that 
with a real singlet, one could either enhance the EWPT or have dark matter, but not both 
simultaneously. One of our main points is that a real singlet is adequate if one relaxes the 
requirement that S be the dominant dark component. We will show that S could constitute 
up to 3% of the total DM mass density, yet still be relevant for direct detection because of its 
large scattering cross section on nuclei. (For other references about the possible EWPT/DM 
connection, see [49]- [62].) 

Because of the complexity of the finite-temperature effective potential that controls the 
dynamics of the EWPT, it is often desirable to do a fully numerical analysis, as for example 
in ref. [63]. But recently it was emphasized [47, 48] that in some situations where the phase 
transition is strong, analytic methods are applicable, which greatly simplify the search for 
working models. This is the case when there is a large barrier already at tree level between 
the electroweak symmetry breaking vacuum {H) = vo/^/2, (S) = 0, and a nearly degenerate 
symmetric one with (H) = 0, (S) = wq. Then the transition can be triggered by rather weak 
thermal corrections, at a temperature significantly lower than the critical vacuum expectation 
value of the Higgs, Vc, which can be relatively close to the zero-temperature one vq. We take 

^Specifically, the DM should be close to half of the Higgs mass in order to have resonant annihilation 
SS — >■ // through Higgs exchange in the S channel, and several large quartic couplings should combine to 
give a small effective coupling A^- 
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advantage of this in the present study, although the solution of transport equations needed 
to compute the baryon asymmetry must still be done numerically. 

In addition to promoting 5 to a dark matter candidate, we take advantage of it to get the 
CP-violation required for baryogenesis by introducing a dimension-6 operator, that modifies 
the top quark mass at nonzero S. The full mass term takes the form 

ytO^i? (l + ^52) f« + h.c. (1.1) 

where r/ is a complex phase and A is a new physics scale. During the EWPT, the top quark 
mass thus gets a spatially-varying complex phase along the bubble wall profile, which provides 
the source of CP violation needed to generate the baryon asymmetry. Ref. [48] considered the 
analogous dimension-5 operator involving S/A, but here we are forced to use S"^ / h? because 
of the Z2 symmetry S — )• —S needed to prevent decay of 5", as befits a dark matter candidate. 

We review the method of construction of the effective potential in section 2, constraints 
from invisible Higgs decays in section 3, and direct detection constraints on the scalar dark 
matter candidate in section 4 along with some results from a random scan over model pa- 
rameters. The absence of other constraints on the model is explained in section 5. The 
computation and resulting distributions of value for the baryon asymmetry are described in 
section 6. Conclusions are given in section 7. 



2 Effective potential 

We follow refs. [47, 48], starting from the tree-level potential for the Higgs doublet H and 
real singlet S*, 

V, = \u{\H\''-\vlf + \K {S'-wlf 

+ \\^\H\''S'' . (2.1) 

This potential has the Z2 symmetry S — t- —S that is needed to guarantee the stability of S as 
a DM particle, but parameters can be chosen such that the Z2 breaks spontaneously at high 
temperatures, giving S a VEV (with H = 0) in the electroweak symmetric vacuum, while the 
true vacuum is along the H axis at T = 0.^ The finite-temperature effective potential for the 
real fields H = /i/\/2 and S can be written in the form 

4 V " u;2 ; 4 
+ 1{T^ -T^){chh^ + CsS^), (2.2) 

where the parameter wq has been traded for its counterpart Wc at the critical temperature of 
the phase transition Tc, Vc is the corresponding critical VEV of /i, and the following relations 
hold: 

vl 

« = Am - 2A/i-|- (2.3) 
wi 

Tl = ^ {vl - vl) . (2.4) 



■^Ref. [48] notes that domain walls associated with this spontaneous breaking of Z2 would only come to 
dominate the energy density of the universe at low temperatures T ~ lO"'^ GeV; but by this time the symmetry 
is restored and the domain walls are no longer present. 
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Here the coefficients Ch and Cs determine the O(T^) corrections to the masses of h and S, and 
are given in terms of the gauge and other couphngs by 

C). = s (99^ + 39" + 129? + \k (2i + "i) + 2") 

=- = A(a.(3|+4|)+2.), (2^5) 

while the zero-temperature masses are given by 

ml = 2\hvl (2.6) 



rric 



hvi + Xh{vl-v'^)(^^-^) . (2.7) 



Counting parameters, and taking the Higgs mass to be determined as m/j = 125 GeV 
[64]- [69], one sees that a given model can be specified by choosing three of them freely, which 
we take to be Am, Vc/wc and vq/vc- It can be shown that as long as vq/vc > 1 so that Tc 
exists, the temperature at which the Z2 spontaneously breaks (so that S gets a VEV prior to 
the EWPT) is always higher than Tc- There is also a restriction 

4 > - (2-8) 

ensuring that the electroweak breaking minimum is the lower of the two at zero temperature. 
More general discussions of the impact of dark matter interacting with the Higgs on stability 
of the electroweak vacuum and perturbativity of the couplings have been given in references 
[39, 71]. In the present work we will enforce perturbativity of at the weak scale by 
restricting the range over which it varies to be less than unity in our Monte Carlo scans. We 
will not make any assumptions about the scale at which the effective theory must be UV- 
completed in order to deal with possible Landau poles, reflecting breakdown of perturbativity 
at higher scales. 

One can efficiently scan the parameter space for models with a strong EWPT charac- 
terized by^ 

^>1, (2.9) 

which is needed to satisfy the sphaleron washout constraint, by randomly varying the inputs 
Am, Vc/wc and vq/v^ over reasonable ranges. This procedure easily generates many examples 
of models with a strong phase transition. However we would like to combine this with the 
requirement that S provides a viable dark matter candidate, which we discuss in the next 
section. 

In addition to the potential (2.2) considered in [47, 48], there is a correction coming from 
the dimension-6 operator (1.1) due to its contribution to the finite-T top quark mass, 

= ^T'^h^{S/Kf (2.10) 
8 

at leading order in the high-T expansion. It does not change the critical temperature as 
defined in (2.2) because it is higher than quadratic order in the fields, and it also does not 

^for caveats concerning the reliability of this estimate, see references [72, 73] 
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change the positions of the critical VEVs since it vanishes, along with its first derivatives, on 
either field axis. The main effect of this operator then is to increase the height of the barrier 
at the critical temperature (making the bubble wall thinner) , and probably to bend the path 
of the bubble wall in field space (to be discussed in section 6). These aspects of the problem 
are being treated in a rather rough way already, and we further require that S/ A remain small 
for consistency, so we do not expect (2.10) to play an important role, and thus omit it from 
our analysis. 

The above procedure implicitly assumes that the high-temperature expansion is valid for 
the models of interest. This is an assumption that can only be checked by redoing the analysis 
using the full effective potential, which could be an interesting undertaking for future work, 
but is beyond the scope of the present paper. We note however that precisely such a check was 
done in ref. [47] to justify this approximation for a general class of singlet Higgs extensions 
to the standard model. There it was found that although the high-T expansion tends to 
somewhat overestimate the strength of the phase transition, nevertheless cases with Vc/Tc > 1 
were always present under the more exact treatment, for models nearby in parameter space 
to those favored by the approximate formalism. 



3 Invisible decay of Higgs boson 

If rris < mfi/2, the S'-particle would contribute to the invisible decay width of the Higgs 
boson. The current constraint on the branching to invisible channels is Br//_^iuvis 

<, 0.35 [70]. 

(The exact value ranges between 0.32 and 0.37 for m/j = 125 — 126.5; our results are not 
sensitive to this variation.) Since the Higgs decay width to visible channels at nih = 125 GeV 
is Fh = 4.07 MeV this imphes 

rj^^5S < 2.2MeV, (3.1) 
while the predicted the decay width for H — t- SS is 

Using (3.2) we can recast (3.1) as 

/ \ 1/4 

A™ ^ 0.036 , '^''^ ) . (3.3) 

\2'^h -rus J 

This constraint is strong enough to effectively exclude all interesting models with rris < mhl'^- 
In the following, we will show that can be no smaller than ~ 0.1 to get a strong phase 
transition. For such values (3.3) then implies \mh — fns < 2 MeV. It will be seen below 
(see for example fig. 2) that there is a narrow, finely-tuned region of parameter space where 
nis — ^mh, corresponding to resonant annihilations of the dark matter in order to achieve low 
enough relic density. This restriction therefore only serves to make this region a bit narrower 
than it would otherwise have been. 



4 Dark matter constraints 

Because of the h'^S'^ cross-coupling, S acquires the interaction X^VohS'^ after electroweak 
symmetry breaking, and so scatters from nuclei by Higgs exchange. Its relic abundance is 
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Figure 1. Comparison of the SS annihilation cross section (in units of the standard rehc density 
value ((Ju)o) as a function of mg using the approximations (4.1) (dashed, blue) and (4.2) (solid, red), 
for the case Am = 0.5. 



partially determined by Higgs-mediated annihilations SS — )• 
there are more important contributions from SS — )• WW, ZZ, 



bb, but if 7715 > ruz^w or rrih, 
hh. Defining rj = mf/m'^, the 



respective contributions to the annihilation cross section from these processes at s 
given by 



\2 



87rm2 [(4 - r/,)2 + rhTl/mj] 
( 6rfe(l - rfe)3/2, 



Ami are 



(4.1) 



X < 



2 + (l-2/r^)2 ^/r^, 



hh 



AA 



rh- 



+ 1 + ^ 



\ 2 

j VI - rh, hh 



where A stands for W or Z, and S^r = 1, 5z = 1/2. However, near nis ~ n-z.w these do not 
include the contributions from 4-body final states due to virtual W and Z emission. A more 
accurate expression for the contributions from all final states except for hh in this region is 
given by factorizing into the SSh fusion part times the virtual h decay using the full width 
of the Higgs, 

2A^r^(s)t;2 



m 



2^2 I ™2p2 



(4.2) 



=4m| 



The hh final state contribution is added to this as in (4.1). In the range 80 GeV< ^/s <340 
GeV we take Tf^{^/s) from ref. [74], which also includes QCD corrections. Above 340 GeV 
we revert to (4.1), because here the 1-loop Higgs self-interaction corrections, also included in 
[74] as for rrih = y/s, would be too large for a light Higgs with nifi = 125 GeV. The resulting 
annihilation cross sections in these approximations are compared in Fig. 1, showing that the 
inclusion of the 4-body final states gives a more smooth dependence on ttt-s near the WW 
threshold. 
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Figure 2. Scatter plot /id vs. ms from a random scan of parameter space with input parameters 
varying in tlie ranges A„i — 0.1 — 1, v^/vc — 1.1 — 10, logj^Q Vc/wc = ~ (+!)• Different groups 
of points are distinguished by their relation to XENON bound: lightest gray points (filled circles) 
are already excluded, orange diamonds ("marginal") are uncertain and blue circles ("allowed") are still 
allowed even by strongest XENON limit. Yellow plus signs show the extension of the allowed region 
when the upper bound on is pushed to 1.5. 



Often an approximation is made by requiring that the total cross section {av) = 
be equal to the standard value {av)^ = 3 x 10~^^cm^/s to get the correct relic density 
from thermal freezeout. We find however, in agreement with [75], that this procedure may 
overestimate relic density by as much as 30-40% in the range of interest. We thus adopt the 
more accurate but still quite efficient freeze-out formalism described for example in refs. [76, 
77]. 

To quantify the DM abundance we define the ratio 

^ 0.11 ^ ' 

In the Monte Carlo scan over models that that selects for a strong EWPT (to be discussed in 
greater detail below), we find that none of them can give /j-ei — 1 while remaining consistent 
with the direct detection constraints (see figure 2). However, there is no need for S to be 
the only DM component; it could make a subdominant contribution to the total DM density, 
while still interacting strongly enough with nuclei to be potentially detectable. If /rei < 1, 
then the relic density of S is suppressed relative to the observed value by /i-gi, and larger 
values of the coupling Am, which controls the barrier height and thus the strength of the 
phase transition, become allowed. We will exploit this possibility in the following. 

Before coming to constraints from direct detection, we point out a subtlety in the DM 
abundance determination. Because the Higgs resonance is very narrow, the cross section at 
the pole is quite large, and the annihlation cross section evaluated at s = 4m| is not an 
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Figure 3. Plot of /id as a function of itls in the case Am = 0.5. Blue dashed hne shows the calculation 
using (crw)s=4m| E^iid red solid line the one using the accurate thermally averaged annhilation cross 
section. 

accurate approximation to the thermally averaged cross section when ms ^ mh/2, where the 
WIMP thermal distribution may be overlapping with the pole. The effect of using the accurate 
thermal annihilation cross section (see ref. [78]) is illustrated in figure 3. In fact neither of 
the results shown there can be trusted in the affected range 50 GeV < nis < fnh/2. Whereas 
{av)g^^^2 is an underestimate, the integrated (av) is an overestimate for the annihilation 
rate, because in reality only a narrow range of momenta see the pole at any given time and 
rapid kinetic-equilibrium-restoring processes would be needed to continuously feed new states 
into this range. We bracketed the uncertainty by performing the computation both ways. In 
the end this does not affect our results because all models sensitive to the effect turn out to 
be excluded by the invisible Higgs decay constraint (3.3), regardless of the way in which the 
relic density was computed. 

Next we consider the constraints from direct detection. The cross section for spin- 
independent scattering on nucleons by virtual Higgs exchange is given by [2] 



where ^ = m„ms/ {mn+nis) is the DM-nucleon reduced mass. The Higgs-nucleon coupling 
suffers from hadronic uncertainties, which have been estimated to be as large as fj\f = 0.26 — 
0.63 [34]. However recent lattice studies obtain roughly consistent values /jv — 0.33 — 0.36 
[79, 80]. For definiteness we adopt the value /tv = 0.35 which also agrees with ref. [81]. 

The XENON 100 experiment currently sets the strongest bounds on agi [82] in the range 
of rus that is of interest for the EWPT. We would like to use their limit to constrain the ratio 
/rei. That limit of course is based upon the assumption /rei = 1, and becomes proportionally 
weaker if /j-ei < 1. We define the fraction 
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Figure 4. Distributions of parameters satisfying the constraints (2.8, 2.9, 3.1) and the nominal 
DM direct detection bound (4.6). Top row shows input parameters, bottom two rows are derived. 
Dimensionful quantities are in GeV units. 



where uxe is the XENONIOO 90% c.l. upper hmit for standard DM of a given mass, and asi 
is the predicted cross section (4.4). Then the constraint on subdominant DM becomes 

/rel < fXe ■ (4.6) 

We find that these bounds can be satisfied for many models that give rise to a strong EWPT. 
For example, a random sample of 2 x 10^ models where the input parameters are varied over 
the ranges \m = 0.1 — 1, vq/vc = 1-1 — 10, log^g Wc/tfc G 1) produces 22500 models 
consistent with the constraint (4.6) as well as with the sphaleron washout bound (2.9), the 
consistency requirement (2.8) and the invisible Higgs decay width (3.1) of previous sections. 
Distributions of various parameters in this set of models can be seen in Fig. 4. One observes 
that the DM mass is typically in the range 80 — 160 GeV, for our choice < 1. (Fig. 2 
illustrates that higher masses are correlated with larger values of A^)- The Vc values fall in 
the range 140 — 220 GeV and as Tc tends to be around 100 GeV strong phase transitions are 
found with Vc/Tc as high as 3.5. The Wc distribution peaks at Wc ~ 160 GeV with Wc < 500 
GeV and the relic density fraction /rei tends to be < 0.01. 

We show the scatter plot of accepted models in /rei versus ms in figure 2 and the same 
data in figure 5 as m-s versus a^s = frcl'^si- The cross section a^s indicates the reach of 



-8- 



T 




Figure 5. Scatter plot of the expected cross section acs = /rci Csj in Xenon experiment vs. 
corresponding to the data shown in figure 2, also with the same color coding. Solid blue ("Nominal 
bound") line shows the nominal XENONIOO bound [82] and the dashed orange lines the more conser- 
vative (upper curve) and more optimistic (lower curve) bounds reflecting uncertainties due to local 
DM distribution. 

the future XENON experiments to rule out a given model, or to verify the existence of its 
associated DM particle. All direct DM bounds inevitably suffer from uncertainties in the 
local Galactic abundance and velocity distribution of the DM. We estimate the effect of 
these uncertainties on the latest XENONIOO constraint following ref. [83], which shows that 
the constraint derived from standard assumptions about the local DM distribution could 
reasonably be stronger or weaker by the respective amounts in Alog^^g'^Xe- 

-0.29128 + 0.4557a; - 0.081349x2, upper 
-1.3329 + 1.4365 x - 0.59337^2 + 0.080595 x^ lower 

where x = log]^o("^s/G'eV). These are our digitizations of the 90% c.l. curves of fig. 5 [left, 
bottom] of ]83]. The altered bounds are shown by the dashed orange lines around the solid 
blue ("Nominal bound") line corresponding to the current XENONIOO bound [82], also shown 
in fig. 5. The colour-coding of the points in Figures 2,5 are the same. The gray regions 
are already excluded by even the weakest of the bounds, the orange regions ("Marginal") fall 
into the region of uncertainty and the blue ones ("Allowed") are still allowed. All of these 
correspond to models with Am < 1, whereas yellow indicates the extension of the allowed 
region when the upper bound on the coupling is relaxed to Xm < 1.5. Note that the current 
XENONIOO results have only begun to encroach on the allowed parameter values. On the 
other hand it is interesting that most of the models (with the exception of the rare ones 
where mg = mh/2) will become discoverable in the relatively near future as direct detection 
sensitivity continues to improve. 
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5 Other constraints 

Because of the inert nature of the S and the fact that it gets a VEV only at high temperature, 
we are free from the additional constraints that had to be imposed on similar models where 
the singlet Higgs was not required to be a DM candidate. For example, nothing prevents 
us from choosing the phase r] in (1.1) to be maximally CP- violating. Ref. [48] considers the 
two-loop Barr-Zee contributions to the electric dipole moments of the electron and neutron. 
But this requires h-s mixing, which does not occur in our model. Ours is similar to models 
in which CP is broken spontaneously at high temperature in this respect. 

Because of the singlet nature of S and its sole couplings being to the Higgs (without 
mixing), and through the dimension-6 operator (1.1), there are no other direct constraints on 
its mass from collider searches, nor from precision electroweak observables. 



6 Baryon asymmetry 

The baryon asymmetry depends upon a source of CP violation that biases sphaleron inter- 
actions near the expanding bubble walls toward production of baryons, as opposed to an- 
tibaryons. We take our relevant CP-violating parameter to be the phase r] in the dimension-6 
coupling in (1.1), and for definiteness we fix 77 = e*'^/^ to maximize the CP violation. (Since 
the baryon asymmetry rjB goes linearly in the imaginary part, the generalization to arbitrary 
phases is straightforward.) Then inside the bubble walls during the phase transition, the top 
quark has a spatially varying complex mass, given by 

mt{z) = ^=h{z) (l + t^^\ = \mt{z)\e''^'^ (6.1) 



V2 ' '\ 

where z is taken to be the coordinate transverse to the wall, in the limit that it has grown large 
enough to be approximated as planar. The existence of the nontrivial phase 9{z) = S{z)/K 
is sufficient to source the baryon asymmetry. In the following, we will initially fix A = 1 TeV 
for the computation r/^. Since r/g ~ 1/A^ for large A, one can always rescale A to adjust 17^ 
to the desired value. 

We follow ref. [48] in approximating the bubble wall profiles in the form 

h{z) = |tic(l + tanh(z/L^)) 

S{z) = \wc{l - tanh(z/L^)) (6.2) 
where the wall thickness is taken to be 

This fully determines the top quark mass profile for a given model. 

The baryon asymmetry is determined by first solving transport equations for the chem- 
ical potentials and velocity perturbations of various fields that develop an asymmetry in the 
vicinity of the bubble wall. We improve upon the treatment given in [48] by using the more 
recent and complete transport equations of [84], which are based on the semiclassical baryo- 
genesis mechanism of refs. [85-89] that determine the chemical potentials of t^, t^, 6^ (the 
left-handed bottom quark) and /i, rather than those of [90]. We also correct an apparent 
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Figure 6. Distributions of wall thickness L^, in units ^, and baryon asymmetry in units 
of the observed value. The latter can also be read as the distribution of (A/TeV)^ for models with 

Vb — ^B.obs- 

error in [48] where there was a mismatch between the orientation of the bubble wall and 
the transport equations that were solved. (The transport equations are not symmetric under 
z — )• — z because it matters whether the wall is expanding into the symmetric phase (correct) 
or into the broken phase (incorrect).) 

One then calculates the left-handed baryon chemical potential 

fiB^ = 5(1 + 4i^i,t)/it^ + 5(1 + ^Kl^b)^^bL + '^Ki^mn > (6-4) 

which is a linear combination of the quark chemical potentials; see Fig. 8 of ref. [63] for for a 
graph of the functions Ki^i{mi{z) /T). The baryon asymmetry depends upon the integral of 
I^Bl o'^^i' the symmetric phase in front of the wall 



405r, 



^ /d^A^B./sphe-^^^-^l^l/^^^™) (6.5) 



where /sph is a function that quickly goes to zero in the broken phase [63], to model the 
spatial dependence of the sphaleron interaction rate, whose value in the symmetric phase is 
Tgph = 10~^ T. We take the wall velocity to be =0.1. Since ~ for small v^, which 
cancels the factor of 1/vw in (6.5), our results are not very sensitive to the value of in the 
expected range ~ 0.01 — 0.1. (For v^j ^ 0.01, the slowly moving wall does not provide 
enough departure from equilibrium to prevent washout by sphaleron interactions, as encoded 
by the exponential factor.) 

We solve the transport equations both by shooting and by relaxation algorithms with 
consistent results. (We find that relaxation gives somewhat more accurate solutions for the 
chemical potential profiles, and we used this method for the results presented here.) Of the 
randomly generated models passing all other constraints, discussed in the previous sections, 
we find that 11% have r]B exceeding the observed value, assuming the scale suppressing the 
dimension-5 source of CP violation is fixed to be A = 1 TeV. This demonstrates that it is 
relatively easy to get sufficient baryogenesis once the phase transition is required to be strong 
enough {vc/Tc > 1). We display the distributions of the wall thickness L^, and 77^ (in units 
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Figure 7. Frequency of wH h? from random scan of models, which should be small to justify the 
effective field theory treatment of the dimension-6 CP-violating operator. 

of the observed value, r/g obs = 8.7 x 10^^^) in Fig. 6. If we rescale A so as to give r/g the 
observed value, then Fig. 6(b) can be reinterpreted as the distribution of values of in TeV^. 

We find that large values of the baryon asymmetry are correlated with strong phase 
transitions, which is expected since the source of CP violation depends partly upon the 
amount by which \mt\ ~ changes in the bubble wall. They also correlate with large values 
of Wc (the critical VEV of S) since this controls the imaginary part of the top quark mass in the 
bubble wall. For consistency however, one should demand that the dimension-6 contribution 
to mt{z) be small compared to the renormalizable one, hence (ifc/A)^ should be small. After 
adjusting A to make r/^ = i]b, ohsi we find that w^/ is distributed as shown in Fig. 7. All 
accepted models have w^j < 1 and 29% of them satisfy w^/ < 0.1. 

From fig. 6 we see that both Ly^ and A tend to be somewhat small. These tendencies are 
actually correlated as is evident from fig. 8. This may be a limitation because for very small 
wall thicknesses our baryogenesis formalism which relies on the semiclassical method [85- 
89] source becomes less reliable. To illustrate the tension, fig. 9 displays the distribution of 
allowed models subjected to additional constraints L^Tc > 3 and A > 500 GeV. Interestingly, 
these cuts tend to remove models with large A^, making the scenario even more testable by 
direct DM searches. But even after these stringent cuts there remains a much greater fraction 
of parameter space with viable examples than was found in the general two-Higgs doublet 
model [63], where a Markov-chain Monte Carlo approach was needed to find a significant 
number of working models. 



7 Conclusions 

In this paper we have established a new possible connection between electroweak baryogenesis 
and dark matter, via a singlet Higgs particle S. Even though S cannot constitute the majority 
of the dark matter in this scenario (at most ~ 3%), it may still be on the verge of discovery 
by XENON-like direct detection searches, due to its large Higgs-mediated cross section on 
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A (GeV) 



Figure 8. Scatter plot A vs. iu,Tc for the models consistent with constraints of sections 2 and 3. 
Dashed lines indicate the consistency constraints A > 500 GeV and > 3T^^. The colour-coding 
of points is the same as in figures 2 and 5. 




Figure 9. Scatter plots corresponding to figures 2 and 5, showing all models that are also consistent 
with the further constraints L^Tc > 3 and A > 500 GeV. 



nucleons. Our predicted range of masses rus = 80 — 160 GeV (for the case where Am is 
restricted to be < 1) is near that where XENONIOO is most sensitive. 

Apart from providing a new dark matter candidate, this theory is interesting in that 
it is easy to find parameters giving a sufficient baryon asymmetry, thanks in part to the 
mechanism of getting a strong electroweak phase transition from having a large potential 
barrier at tree level, that allows an analytic treatment of the phase transition properties. The 
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model is not UV-complete because we have invoked a nonrenormalizable operator Q^HS^tji 
to get CP violation. However such an operator could be obtained from a renormalizable 
theory by integrating out a heavy Higgs doublet H' that does not get a VEV, has a Yukawa 
coupling to the top quark and a (possibly complex) H^H'S^ quartic coupling. We found 
that the new physics scale A could be sufficiently large to justify the low-energy effective field 
theory treatment. 

It would be interesting to repeat this exercise in the context of two-Higgs doublet models, 
where the second doublet is not very heavy. In the purely numerical search done in ref. [63], 
there is not much insight into what is special about the models that work (other than that 
large values of S are needed in the bubble wall to get enough baryon asymmetry, in common 
with the present situation). It is possible that they are variants of the mechanism used here, 
involving the large tree-level potential barrier. Moreover, the studies of the Inert Doublet 
Model done in ref. [17-19] could likely be generalized to find a similar DM mass range as we 
found here, if the possibility of subdominant DM and all the relevant annihilation channels 
were included. 

It is intriguing that the DM mass range we found includes 130 GeV, since there are 
exciting indications of DM annihilation into two photons of approximately this energy [92]- 
[96] . Recently it was proposed that subdominant DM could account for this observation [97] . 
In that case it was necessary to have both a larger fraction of DM in the 130 GeV component, 
and for the dominant annihilation channel to be into two photons, neither of which are the 
case in the present scenario. However these conclusions may not necessarily hold for other 
models, and so it could be worthwhile to explore whether the phenomena we focus on here 
could be compatible with a simultaneous explanation of the 130 GeV line. 
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